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Abstract. The first ever double pulsar, discovered by our team a few 
months ago, consists of two pulsars, one with period of 22 ms and the 
other with a period of 2.7 s. This binary system with a period of only 2.4- 
hr provides a truly unique laboratory for relativistic gravitational physics. 
In this contribution we summarize the published results and look at the 
prospects of future observations. 



1. Introduction 

In these proceedings, the contribution by Burgay et al. describes the discovery 
of the compact double neutron star system J0737— 3039 and its implication for 
the detection rate for ground-based gravitational wave detectors and tests of 
general relativity (see also Burgay et al. 2003). With the discovery of radio 
pulses from the companion, the system becomes unique in its potential as a 
testbed for relativistic gravity. The data resulting in the timing parameters 
presented for both pulsars by Lyne et al. (2004) and shown in Table 1, were 
obtained with the Parkes radio telescope and the Lovell telescope at Jodrell 
Bank. An updated timing solution covering a full year of observations will be 
presented in a forthcoming paper (Kramer et al. in prep.). The contribution by 
Manchester et al. discusses the implications for studying pulsar magnetospheres. 



2. Tests of general relativity 

Since neutron stars are very compact massive objects, double neutron star (DNS) 
binaries can be considered as almost ideal point sources for testing theories of 
gravity in the strong-gravitational-field limit. Tests can be performed when a 
number of relativistic corrections to the Keplerian description of the orbit, the 
so-called "post-Keplerian" (PK) parameters, can be measured. For point masses 
with negligible spin contributions, the PK parameters in each theory should only 
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be functions of the a priori unknown neutron star masses, and rag, and the 
well measurable Keplerian parameters. With the two masses as the only free 
parameters, the measurement of three or more PK parameters over-constrains 
the system, and thereby provides a test ground for theories of gravity. In a theory 
that describes the binary system correctly, the PK parameters produce theory- 
dependent lines in a mass-mass diagram that all intersect in a single point. (See 
contributions by Weisberg & Taylor for an update on the tests provided by PSR 
B1913+16 and by Stairs et al. for tests possible with PSR B1534+12.) As shown 
in Table 1, we have measured PSR J0737— 3037A's advance of periastron, u>, and 
the gravitational redshift/time dilation parameter 7 and have also detected the 
Shapiro delay in the pulse arrival times of A due to the gravitational field of 
PSR J0737— 3039B (see Fig. 1), providing a precise measurement of the orbital 
inclination of s = s'mi = 0.9995(lg 2 ) and the "range" parameter r. This already 
provides four measured PK parameters, resulting in a -mg plot shown in 
Figure 2 where we label cb(mA,mB), 7^, tub), r(m J 4,ms), s(mA,mB) as 
predicted by general relativity (Damour & Deruelle 1986). 

In addition to tests with these PK parameters, the detection of B as a pulsar 
opens up opportunities that go well beyond what is possible with previously 
known DNSs. With a measurement of the projected semi-major axes of the 
orbits of both A and B (see Table 1), we obtain a precise measurement of the 
mass ratio, R(m\, me) = m\/mB = xb/xa, providing a further constraint 
displayed in Figure 2. For every realistic theory of gravity, we can expect the 
mass ratio, R, to follow this simple relation (Damour & Taylor 1992). Most 
importantly, the i?-line is not only theory- independent, but also independent 
of strong- field (self-field) effects which is not the case for PK-parameters. This 
provides a stringent and new constraint for tests of gravitational theories as any 
intersection of the PK-parameters must be located on the i?-line in Fig. 2. With 
four PK parameters already available, this additional constraint makes J0737- 
3039 the most overdetermined DNS system to date where the most relativistic 
effects can be studied in the strong-field limit. 

3. Future potential 

We expect to measure the orbital period derivative, P^, by autumn 2004. We 
note that the distance to this pulsar will be needed to correctly account for 
relative accelerations and hence to interpret i\ in the context of gravitational 
theories. VLBI observations to determine the system's parallax are underway. 
In the longer term, further effects should be measurable. 

3.1. Geodetic Precession & Aberration 

In general relativity, the proper reference frame of a freely-falling object suffers 
a precession with respect to a distant observer, called geodetic precession. In 
a binary pulsar system this geodetic precession leads to a relativistic spin-orbit 
coupling, analogous of spin-orbit coupling in atomic physics (Damour & Ruffini 
1974). As a consequence, the pulsar spin precesses about the total angular 
momentum, changing the relative orientation of the pulsar towards Earth. Since 
the orbital angular momentum is much larger than the pulsars' spins, the total 
angular momentum is practically represented by the orbital spin. The precession 
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Table 1. Observed and derived parameters of PSRs J0737-3039A 
and B. Standard errors are given in parentheses after the values and 
are in units of the least significant digit (s). 



Pulsar 

Pulse period P (ms) 
Period derivative P 
Epoch of period (MJD) 
Right ascension (J2000) 
Declination (J2000) 
Orbital period Ph (day) 
Eccentricity e 

Epoch of periastron To (MJD) 

Longitude of periastron u> (deg) 

Projected semi-major axis x = a sin i/c (sec) 

Advance of periastron w (deg/yr) 

Gravitational redshift parameter 7 (ms) 

Shapiro delay parameter s 

Shapiro delay parameter r (/is) 

Total system mass m\ + rriB (Mq) 

Mass ratio R = m^/m-B 

Orbital inclination from Shapiro s (deg) 

Orbital inclination from (R, u>) (deg) 

Stellar mass from (R, ui) (Mq) 



PSR J0737-3039A PSR J0737-3039B 
22.69937855615(6) 2773.4607474(4) 
1.74(5) x 1(T 1S 0.88(13) x 10~ 15 
52870.0 52870.0 
07 h 37 m 51 s . 247(2) 
-30°39'40".74(3) 
0.102251563(1) 
0.087779(5) 
52870.0120589(6) 

73.805(3) 73.805 + 180.0 

1.41504(2) 1.513(4) 
16.90(1) 
0.38(5) 
0.9995(-32, +4) 
5.6(-12, +18) 

2.588(3) 
1.069(6) 
87(3) 
87.7(-29, +17) 
1.337(5) 1.250(5) 



rate (Barker & O'Connell 1975, Borner et al. 1975) depends on the period and 
the eccentricity of the orbit as well as the pulsar and companion mass. With 
the parameters shown in Table 1, general relativity predicts precession periods 
of only 75 yr for A and 71 yr for B. 

Geodetic precession should have a direct effect on the timing as it causes 
the polar angles of the spins and hence the effects of aberration to change with 
time (Damour & Taylor 1992). These changes modify the observed orbital pa- 
rameters, like projected semi-major axis and eccentricity, which differ from the 
intrinsic values by an aberration dependent term, potentially allowing us to 
infer the system geometry. Other consequences of geodetic precession can be 
expected to be detected much sooner. These arise from changes in the pulse 
shape and its polarisation properties due to changing cuts through the emission 
beam as the pulsar spin axis precesses. Due to these effects, which may compli- 
cate the TOA determination, geodetic precession is detected for PSRs B1913+16 
(Kramer 1998; Weisberg & Taylor, these proceedings) and B1534+12 (Arzou- 
manian 1995; Stairs et al., these proceedings). For PSR B1913+16 it is used to 
derive a two-dimensional map of the emission beam (Weisberg & Taylor 2002) 
and to determine the full system geometry leading to the prediction that the 
pulsar will not be observable from Earth around 2025 (Kramer 1998). Since the 
precession rates for A and B are about four times larger, the visible effects in 
J0737— 3039 should be much easier to detect. 
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Figure 1 . The effect of the Shapiro delay caused by the gravitational 
potential of B seen in the timing residuals of A. Residuals were aver- 
aged into 75 equal bins of orbital phase. Top: timing residuals obtained 
by fitting all model parameters shown in Table 1 except the Shapiro 
delay parameters r and s. The left-over structure represents the higher 
harmonics of the Shapiro delay that are unabsorbed by fits to the Roe- 
mer delay. Middle: Timing residuals after fitting additionally for the 
Shapiro delay. Bottom: as middle display, but with the Shapiro delay 
parameters r and s set to zero. 




3.2. Higher Post-Newtonian orders 

In all previous tests of general relativity involving PK parameters, it was suf- 
ficient to compare the observed values to theoretical expressions computed to 
their lowest post-Newtonian approximation. However, higher-order corrections 
may become important if relativistic effects are large and timing precision is 
sufficiently high. Whilst this has not been the case in the past, the double pul- 
sar system may allow measurements of these effects in the future (see Lyne et 
al. 2004). 

Light propagation effects Higher-order corrections to the present description 
of the light propagation would account for the assumptions currently made in 
the computation of the Shapiro delay that gravitational potentials are static 
and weak everywhere. Corresponding calculations taking these next higher or- 
der corrections into account have been presented by Wex (1995a), Kopeikin & 
Schafer (1999) and Kopeikin (2003). 

In addition to light propagation effects causing time delays, others exist that 
are related to light bending and its consequences. Light bending as discussed by 
Doroshenko & Kopeikin (1995) would be superposed on the Shapiro delay as a 
typically much weaker signal, which arises due to a modulation of the pulsars' 
rotational phase by the effect of gravitational deflection of the light in the field 
of the pulsar's companion. This aberration-like effect would only be visible for 
very nearly edge-on systems since it depends crucially on the orientation of the 
pulsar's spin-axis in space. The double pulsar system fulfils these requirements. 
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Figure 2. The observational constraints upon the masses toa and 
niB- The solid regions are those which are excluded by the Keple- 
rian mass functions of the two pulsars. Four further constraints are 
shown: (a) the measurement of the advance of periastron to; (b) the 
measurement of R = me/mA; (c) gravitational redshift parameter 7 
and (d) Shapiro delay parameters r and s. the intersection of the four 
constraints, with the scales increased by a factor of 50. The permitted 
regions are those between the pairs of parallel lines and the solid region 
is the only area permitted by all constraints. 



Moment of inertia In contrast to Newtonian physics, general relativity predicts 
that the neutron stars' spins affect their orbital motion via spin-orbit coupling. 
This effect would be most clearly visible as a contribution to the observed Cj in a 
secular (Barker & O'Connell 1975), and periodic fashion (Wex 1995b). For the 
J0737— 3039 system, the expected contribution is about an order of magnitude 
larger than for PSR B1913+16 (Lyne et al. 2004), i.e. 2 x 10~ 4 deg yr^ 1 (for A, 
assuming a geometry as determined for PSR B1913+16; Kramer 1998). As the 
exact value depends on the pulsars' moment of inertia, a potential measurement 
of this effect would allow us the moment of inertia of a neutron star to be 
determined for the first time (Damour & Schafer 1988). If two parameters, 
e.g. the Shapiro parameter s and the mass ratio R, can be measured sufficiently 
accurately, an expected Co eyLp can be computed from the intersection point. This 
value can be compared to the observed value w b s which is given by (see Damour 
& Schafer 1988) 



1 + Aw 2PN - /A^) - 9 B ^Eo\ , (1) 

where the last two terms represent contributions from the pulsars' spins. In these 
terms, g A ' B are geometry-dependent factors whilst Aw^g B arise from relativistic 
spin-orbit coupling (Barker & O'Connell 1975), formally at the level of first post- 
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Newtonian (1PN) approximation. However, it turns out that for binary pulsars 
these effects have a magnitude equivalent to second post-Newtonian (2PN) order 
(Wex 1995b), so that they only need to be considered if Co is to be studied at 
this higher level of approximation. We find Awso oc I/Pm 2 (Damour & Schafer 
1988), so that with precisely measured masses m the moment of inertia I can 
be measured and the neutron star equation-of-state and our understanding of 
matter at extreme pressure and densities can be tested. 

The dependence of Awgo on the spin period P suggests that only a mea- 
surement for pulsar A can be obtained. It also requires that at least two other 
parameters can be measured to a similar accuracy as Co. Despite being a tough 
challenge, e.g. due to the expected profile variation caused by geodetic preces- 
sion, the prospects are promising. Simulations indicate that a few years of high 
precision timing are sufficient (Kramer et al. in prep.). 

4. Conclusions 

The double pulsar system offers improved but also new tests of general relativ- 
ity. The current data indicate an agreement of the observed with the expected 
Shapiro parameter of s bs/s CX p = l.OOOlzLO. 00220 (Kramer et al. in prep.) where 
the uncertainties are likely to decrease, enabling the other discussed applications. 

Acknowledgments. MK thanks the Aspen Center for Physics and the or- 
ganisers for the received financial support. 
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